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ABSTRACT 

Equations of light, propagating from quasar to observer on earth, are integrated in 
the time-dependent gravitational field of the solar system by making use of either 
retarded or advanced solutions of the Einstein field equations. This technique 
allows to separate explicitly the effects associated with the propagation of gravity 
from those associated with light in the integral expression for the relativistic time 
delay of light. We prove that the relativistic correction to the Shapiro time delay 
discovered by Kopeikin (ApJ, 556, LI, 2001) changes sign if one retains direction 
of the light propagation but replaces retarded with advanced solution of the 
gravitational field equations. Hence, this correction is due to the propagation of 
gravity and allows to measure its speed. Relativistic VLBI experiment conducted 
by Fomalont and Kopeikin in September 2002 proves that the propagation of 
gravitational field is characterized by the retarded potentials. 

Subject headings: gravitation - relativity - techniques: interferometric - galaxies: 
quasars 



1. Introduction 

Einstein's equations of general relativity have a hyperbolic character (Landau & Lifshitz 
1971; Wald 1984). From mathematical point of view a hyperbolic wave equation allows in 
general two types of solutions described respectively by retarded and advanced potentials 
(Synge 1960; Aharoni 1985). It is often supposed that the advanced solution is unphysical 
because it violates the physical principle of causality. However, the theory of general relativ- 
ity does not exclude the advanced solution per se. Making a choice between the retarded and 
advanced solutions is a matter of a boundary condition imposed on the Einstein equations. 
Although the causality-preserving boundary condition is reasonable from theoretical point 
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of view, one does not know a priori what kind of condition is realized for the gravitational 
field in real world. For this reason, experimental study of the nature of the boundary condi- 
tion for Einstein's equations has a great importance for foundations of theoretical physics. 
Timing observations of binary pulsar PSR 1913+16 indicates that the retarded solution of 
the Einstein equations is the most appropriate for the description of orbital dynamics of 
the neutron stars (Weisberg & Taylor 2002). Nevertheless, an independent experimental 
verification of the retarded nature of the gravity propagation is required. 

In this paper we solve null geodesic equations for hght propagating from a distant source 
(quasar) in time-dependent gravitational field of moving bodies (solar system) by making 
use of the mathematical technique invented by Kopeikin (1997) and further elaborated in 
(Kopeikin & Schafer 1999; Kopeikin & Mashhoon 2002). We derive the relativistic delay in 
the time of propagation of light and analyze how it relates to the retarded and advanced 
solutions of the Einstein equations. This give us an unambiguous criterion for distinguishing 
the relativistic effects associated with the propagation of gravity from those related to the 
propagation of light - the question which was a matter of confusion for some researchers 
(Asada 2002; Will 2003). 

In this paper we assume that the propagation of light from quasar towards observer 
always goes along the electromagnetic null-cone from the past to the future. On the other 
hand, we allow the wave solutions of the Einstein equations describing propagation of grav- 
itational field to be both retarded and advanced. Difference between them is incorporated 
to the time variable that (for each solution) is an implicit function of the coordinate time 
t given by the gravity null-cone equation 



where a; is a field point, Xa{s^) is position of ath gravitating body ^ (a = 1, 2, A^) depending 
on the time argument s^- Numerical parameter e = —c/cg, where c is the limiting speed of 
special relativity and Cg is the 'speed-of-gravity' parameter taking either positive or negative 
values such that in general relativity \cg\ = c. 

In order to distinguish the gravity-propagation effects from those associated with light 
we shall operate separately with the advanced gravity null cone characterized formally by 
the condition —Cg = c, and with the retarded gravity null cone singled out by the condition 
Cg = c. Thus, the parameter e in Eq. (1) takes two different numerical values 
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retarded time , 
advanced time , 
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Notice that the variable is different for each body as it depends on the body's position Xa as weh. 
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such that = 1. The parameter e plays a supplementary role helping to track down the 
presence of the gravity propagation effects in the solutions of both Einstein's and the light- 
ray-propagation equations. This tracking property of the parameter e is used for providing 
unique and unambiguous physical interpretation of the relativistic effect in the VLBI time 
delay predicted by Kopeikin (2001) and observed by Fomalont & Kopeikin (2003). 

In what follows we rely upon the system of units such that G = c = 1. Boldface italic 
letters denote spatial vectors, e.g, a = = {a^,a^, a^). The dot between two spatial vectors 
is the Euclidean scalar product, a-b = a^b^ + a^b^ + a^b^. The cross between two spatial 
vectors denotes the Euclidean vector product, ax b = Sijka^b'^, where Eijk is the Levi-Civita 
symbol such that £123 = +1. Roman indices run from 1 to 3 while Greek indices run from 
to 3. Repeated indices assume the Einstein summation rule. Greek indices are raised and 
lowered with the Minkowski metric rjap = diag(— 1, +1, +1, +1). 



2. The Retarded and Advanced Solutions of the Einstein Equations 

We assume that the gravitational field is weak everywhere so that the metric tensor Qajsit, x) = 
Vai3 + hai3(t,x), where hap{t,x) is the perturbation of the gravitational field, and t and x 
are time and spatial coordinates. 

Let gravitational field be generated by the stress-energy tensor of moving point-like 
massive particles (Landau & Lifshitz 1971) 

N 

T^^{t,x) = Y.M,X\t)u:{t)u^,{t)6{x-Xa{t)), (3) 

a=l 

where a = 1, 2, enumerates gravitating bodies of the solar system. Ma is the (constant) 
rest mass of the ath body, Xa(t) are time-dependent spatial coordinates of the ath body, 
Vait) = dxa{t)/dt is the orbital velocity of the ath body, m° = IIa(l, Va) is the four-velocity 
of the ath body, = (1 — f^)"^''^, and S{x) is a 3-dimensional Dirac's delta-function. 

We shall work in the harmonic coordinates (Weinberg 1972) with the origin at the 
barycenter of the solar system. Formal mathematical solution of the Einstein equations in 
the harmonic coordinates and in the linearized post-Minkowskian approximation (Bel at al. 
1981; Damour 1987) is given by the Lienard-Wiechert tensor potential 

where ra{se) = x — Xa{se), ra{se) = \ra{se)\, Va{se) = dxa{se)/dse, and the time argument 
= s^{t,x) is determined as an implicit solution of the retarded (e = —1) or advanced 
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(e = +1) gravity null-cone Eq. (1). At this stage it is obvious that because we were solving 
the Einstein equations only gravity propagation with the speed c is involved in Eq. (4) 
through the retarded (advanced) time argument s^. 

Non-linear corrections to the linear solution (4) of the Einstein equations can be obtained 
by making use of the iterative post-Minkowskian approximations (Bel at al. 1981; Damour 
1987). This non-linear metric tensor should be substituted to the equations of motion of the 
bodies comprising the gravitating system under consideration. Solution of the equations of 
motion gives us the parametric trajectory Xa{t) of each body orbiting around the barycenter 
of the iV-body system with the instantaneous velocity Va{t) = dxa{t) / dt, and the acceleration 
aa{t) = dva{t)/dt. We shall assume that such solution of the body's equations of motion 
has been found and used in the expression for the metric tensor perturbation given by Eq. 
(4). More involved discussion of this problem in the framework of the post-Minkowskian 
approximations is given by Bel at al. (1981) and Damour (1987). 



3. Light-Ray Trajectory Parameterization 

We consider motion of light in the non- stationary gravitational field described by the metric 
Qajj and assume that light has no back action on the gravitational field. Hence, one can 
use equations of hght geodesies. Let the motion of light be defined by the initial-boundary 
conditions 

, , dxi—oo) , , , 

xih) = xo , — = k , (5) 

where fc^ = 1. These conditions define the coordinates Xo of the photon at the instant of 
emission of light, to, and direction of its velocity at the past null infinity (Wald 1984). 

The original null geodesic equations of light propagation are rather complicated. These 
equations can be simplified and recast to the integrable form (Kopeikin & Schafer 1999; 
Kopeikin & Mashhoon 2002). In the case of weak gravitational field the integration is 
performed by making use of iterations (Kopeikin 1997). One starts from the unperturbed 
trajectory of the light ray that is a straight line 

x'{t) = x%{t)=xi + k'{t-to) , (6) 

where the coordinate time t is a running parameter along the light ray, and tg, Xq, = k 
have been defined in Eq. (5). In this approximation velocity of the photon is = k^, and is 
considered as constant in the next iteration involving the metric tensor perturbation hap- 

It is convenient to introduce a new independent parameter a along the photon's trajec- 
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tory according to the rule 

a = k ■ x^it) = t — to + k ■ xq . (7) 

The time to of the hght emission corresponds to the numerical value of the parameter ctq = 
k ■ Xq, and the numerical value of the parameter a = corresponds to the instant 

r = to-k-xo, (8) 

that is the time of the closest approach of the unperturbed trajectory of light ray to the 
origin of the coordinate system. We emphasize that the numerical value of the moment t* is 
constant for a chosen trajectory of light ray and depends only on the space-time coordinates 
of the point of emission of the photon and direction of its propagation. Thus, we find the 
relationship 

t = f + a, (9) 

which reveals that the variable a is negative from the point of emission up to the point of 
the closest approach, x^{t*) = to the origin of the coordinate system, and is positive 
otherwise. The differential identity dt = da is valid and, for this reason, the integration 
along the light ray's path with respect to time t can be replaced by the integration with 
respect to variable a. 

Making use of the parameter a, the equation of the unperturbed trajectory of the light 
ray can be represented as 

x\a) = x%{a) = k'a + C, (10) 

and the distance, r{a) = |a;7v(^)|; of the photon from the origin of the coordinate system 
reads 

r(a) = V(r^ + d^ . (11) 

The constant vector = $, = kx (xq x k) = k x [xN{t) x k) is called the impact parameter 
of the unperturbed trajectory of the light ray with respect to the origin of the coordinate 
system, d = |^| is the (Euclidean) length of the impact parameter. We notice that the vector 
^ is transverse to the vector k in the Euclidean sense, that is k-$, = 0, where "•" denotes the 
Euclidean dot product. Vector is directed from the origin of the coordinate system towards 
the point of the closest approach of the unperturbed path of the light ray to this origin. This 
vector plays auxiliary role in our discussion and has no essential physical meaning as it can 
be changed by the shift of the origin of the coordinates. Time of the closest approach t* is 
a convenient parameter being useful in the following calculations but it does not enter the 
final formula for the relativistic time delay as it should be expected on the general physical 
ground. 
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4. Relativistic Time Delay 

Time of propagation of light from the point Xq to x\ = x*(ti) is (Kopeikin 2003) 

ti-to = \xi- xo\ + A{ti,to) , (12) 

where ^ 

A(ti, to) = I k'^k" I ' Kp (t, XN{t)) dt , (13) 

^ Jto 

ti is the time of observation, and the (Minkowskian) null vector of the light propagation 
k" = (1,A;*) is directed from the past to the future. For the sake of simplicity we shall 
work out the subsequent calculations only up to the linear terms with respect to the (non- 
constant) velocity Va of the moving bodies. This is presently enough for the purposes of the 
relativistic astrometry in the solar system. 

Substitution of Eq. (4) to (13) and replacement of the variable t ^ a, yields 

where ^^(se) = ^ + ka — a;a(se), and the variable relates to the coordinate time a via Eq. 
(1) of the gravity null cone taken on the unperturbed light ray trajectory (electromagnetic 
null cone) at the point at which the light particle is at the time a 

a + t* = s, - + ka - Xais,)\ . (15) 

We remind once again that the parameter e = ±1 indicates the presence of the gravity- 
propagation effects. 

Differentiation of Eq. (15) yields an exact relationship between the differentials of the 
time variables a and 

da ds, 

= ; 5 (lo) 

ra + eva -Va Ta + ek- r„ 

where = xi^{t) — Xa{s^, and the denominator in the left side of this equation coincides 
exactly with the denominator of the metric tensor perturbation in Eq. (4). Transformation 
(16) allows to replace integration with respect to coordinate time a to that with respect to 
the retarded (advanced) time which is simpler and more informative. 

In order to take integral (14) one introduces a new variable y = Ta + ek-Va and calculates 
it on the gravity null cone Eq. (15). 

ye = ra + ek-ra = -e {t* + k ■ Xa{s,) - s^) ■ (17) 
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We get 

dy^ = e{l - k-Va{s^))ds^ 
so that the above integral (14) is recast to 



N 



a=l 



a=i -fy^o 

Integration by parts yields 

A(ti,to) = 2eJ]M, (l-A;-'u„(s,i))lny,i- (l-fc-i;,(s,o))ln 



N 



a=l 



where 



VeO 



Via + efc ■ Via 

+ ek ■ roa 



(19) 



(20) 



(21) 
(22) 



with ria = xi — a3a(sei), = xq — Xa{seo), and the retarded (advanced) times s^i and s^o 
are connected to the time of observation, ti, and emission of hght, to, via the gravity null 
cone Eq. (1). Eq. (20) does not assume that the velocities of the gravitating bodies are 
constant. It only means that we have neglected accelerations when performing the integral. 



5. Differential VLBI time delay 

Let us assume now that there are two earth-based VLBI stations and the front of an electro- 
magnetic wave from quasar propagates towards the earth. Quasars are almost at the edge of 
the visible universe and they do not reveal any annual parallax and/or proper motion (Wal- 
ter & Sovers 2000). For the array of the VLBI stations spread out across the earth one can 
assume that the wave front of quasar's electromagnetic radiation is plane and characterized 
by a single wave vector k°' = {l,k). Taking two rays from the wave front, directed from the 
quasar towards the first and second VLBI stations, and subtracting Eq. (20) for the first 
light ray from that for the second ray yields 

t2-ti = -K-B + A{t,,t2), (23) 

where A(ti,t2) = A(t2,io) — A(ti,to), K = —k is the unit vector from the barycenter of 
the solar system to the quasar (quasar's annual parallax is neglected), and B = X2{ti) — 
Xi{ti) is the barycentric baseline between the two VLBI stations. Here ti and t2 are the 
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barycentric coordinate times of arrival of the electromagnetic signal from the quasar to the 
first and second VLBI stations respectively, and Xi and X2 are spatial coordinates of the 
first and second VLBI stations with respect to the barycentric frame of the solar system 
that is chosen as the primary non-rotating reference frame (Walter & Sovers 2000). We 
omitted in Eq. (23) the terms which are proportional to the velocity of motion of the second 
station as these terms are not crucial in discussion pertained to the interpretation of the 
relativistic VLBI experiment conducted on September 8, 2002 (Fomalont & Kopeikin 2003). 
Further theoretical details of the relativistic algorithm of transformation of the barycentric 
to geocentric coordinates are contained in (McCarthy 2000; Walter & Sovers 2000). They 
are not so important in the discussion that follows, but are essential for understanding how 
VLBI works in general. 

The difference A{ti,t2) is obtained from Eqs. (20), (23) after long and tedious calcula- 
tions accompanied by analysis of the residual terms similar to that done in (Kopeikin 2001). 
We have proved that with accuracy better than 1 ps. 



where Va{sei) is the velocity of the ath gravitating body at the retarded (advanced) time 

Sel, = \ria\, = \r2a\, r'la(s,l) = Xi{ti) - £C„(s,i) and r2a{Se2) = X2{t2) - Xa{s,2) ] 

moreover, the retarded (advanced) times s^i and s^2 are calculated according to Eq. (1) 



where again we notice the presence of the gravity-propagation parameter e = ±1. Thus, 
Eqs. (25), (26) reveal the true physical reason for the coordinates Xa of the bodies in Eq. 
(24) to be displaced from their present positions - the finite speed of propagation of gravity. 

Measuring relativistic time delay of radio signal propagating from quasar through the 
gravitational field of moving body allows to evaluate position of the body in the sky at the 
time of observation and compare it with that one would observe if position of the body were 
measured directly in optics by receiving light emitted by the body itself (Seidelmann 1992). 
In the case of the retarded solution of the Einstein equations (e = —1) position of the body 
measured directly in optics and that measured indirectly through the realtivistic time delay 
coincide precisely, that is both optical and gravitational 'images' of the body experience the 
same aberration in complete agreement with general relativity. It is worth noting that for 
slowly-moving particles this effect can not be observed in the approximation under discussion 
because of a different structure of the perturbing gravitational force leading to the apparent 




(24) 



Sel = ti + e|a?i(ti) - Xa{s^i)\ , 

Se2 = t2 + e\X2it2) - Xais,2)\ , 



(25) 
(26) 
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cancellation of retarded effects in the solutions of the post-Newtonian equations of motion of 
the particles (Carlip 2000). This cancellation does not happen for light and the retardation 
of gravity can be observed experimentally by making use of high-precision VLBI technique 
(Fomalont & Kopeikin 2003). 

6. Retarded versus Advanced Solutions of the Einstein Equations in the VLBI 



Let us elaborate Eq. (24) for two different cases corresponding to the retarded (e = — 1) 
and advanced (e = +1) solutions of the Einstein equations. This will help us to split the 
gravity-propagation effects from those associate with the propagation of light. The gravity 
propagation effects depend on the parameter e and changes sign when one uses retarded 
solution instead of advanced one and vice versa. 



The retarded case of e = — 1 has been already treated in the papers (Kopeikin 2001, 2003) 
from various points of view. The effect of retardation of gravity appears in Eq. (24) as a 
displacement of the light-ray deflecting bodies from their present Xa(ti) to retarded positions 
Xa{s_i) {i = 1,2). We notice that velocities of the gravitating bodies, Va, are small with 
respect to the fundamental constant c. Furthermore, the time taken by light ray to cross 
the solar system is much smaller than orbital periods of Sun and other planets around the 
barycenter of the solar system. Thus, one can expand the retarded positions Xa{s^i) of the 
bodies in the Taylor series around their present positions Xa{ti) taken at the arrival times 
ti to the first or second VLBI stations. Moreover, we shall assume that only the impact 
parameter dj of the light ray with respect to Jupiter is small. Hence (Kopeikin 2001), 



where hereafter the subscript J refers to Jupiter, n is the unit vector of the impact parameter 
of the light ray with respect to Jupiter, and 6' is a small angle in the sky between the directions 
towards the undisturbed geometric position of the quasar and that of Jupiter. 

After expanding body's retarded position Xa{s-i) around the instant ti, neglecting all 
terms quadratic with respect to velocities of the gravitating bodies, proportional to their 



Time Delay 



6.1. Retarded Solution 




(27) 
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accelerations and the products of body's velocities with the angle 6, one gets 

a=l ^ ' 

where we have restored the constants G and c for the sake of convenience, Ria = Xi{ti)—Xa(ti) 
is the difference between coordinates of the ith VLBI station and ath gravitating body taken 
at the time of arrival of radio signal to the zth station, Va = Va{ti), and all quantities in the 
second term on the right side of Eq. (28) are also evaluated at time ti. 

The second term in Eq. (28) was obtained as a consequence of expansion of the retarded 
time s_ (e = —1) in the gravity null cone Eqs. (25), (26) which describe for this value of e the 
effect of retardation of gravity. We emphasize that light from the quasar does not propagate 
along the gravity null cone that connects the body (Jupiter, Sun, etc.) and the light particle. 
Light propagates along the null world line, connecting the quasar and observer, which is not 
a part of the gravity null cone ^. Consequently, the origin of the second term in Eq. (28) is 
due to the propagation of gravity as it was shown by Kopeikin (2001). 

To confirm our interpretation from different prospect we shall study the differential 
VLBI time delay for the case of the advanced solution of the Einstein equations. If the 
second term in the right side of Eq. (28) is due to the propagation of gravity, it must change 
sign in the advanced solution. This indeed takes place. 



6.2. Advanced Solution 

Let us now take the advanced solution of Einstein equations, that is e = +1, and calculate 
the impact of the advanced gravitational field on the VLBI time delay. We shall first have 
from Eq. (24) 

A(ti,t,) = -2f2Ma(l + K-v.{s^,)) 1, ^1^(^+0 -^-^i^l^+i) ^ (29) 

^ ^ ^ r2a{s+2) - K-r2a{s+2) 

where the advanced times s+i and s+2 are computed from Eqs. (25), (26) with e 
notice that the logarithmic function in Eq. (29) can be recast as follows 

_ ria(g+i) - K-ria{s+i) _ ria(s+i) + K-ria{s+i) _ 2 |K x ria(s+i)| 

r2a{s+2) - K-r2a{s+2) r2a(s+2) + -f^-T'2a(s+2) \K X r2a{s+2)\ 



^This might be a reason for the misconceptions presented in (Asada 2002; Will 2003) and discussed by 
Kopeikin (2003) in more detail. 
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By making use of the Taylor expansion of Xa{s+2) around the instant s+i one can get 



\Kxr2a{s+2)\ |K X ri„(s+i)|2 

where 

S+2-S+1 = -K-B + N^a-B + 0{va) . (32) 

The second term in the right side of Eq. (31) is neghgibly small for any body of the 
solar system but might be large enough to be measured in the field of Jupiter. For Jupiter, 
one has 

|Kxru(s+i)| _ A{K-B){vj.n) 
''^Kxr,,(^+2)| " r,je ' ^^^^ 

where we put together two terms in Eq. (32) after making use of Eq. (27). One can see 
from Eq. (33) that the second term in the right side of Eq. (30) is by the factor 9 smaller 
than the second term in the right side of Eq. (28). During the time of the VLBI experiment 
(Fomalont & Kopeikin 2003) 6' ~ 3.7 arcminutes and the second term in the right side of Eq. 
(28) was about 51 /ias. Hence, the magnitude of the second term in the right side of Eq.(30) 
is much less than 1 yuas and can be neglected because the precision of the VLBI experiment 
on September 8, 2002 was ~ 10 /las (Fomalont & Kopeikin 2003). 

Thus, we obtain that in the case of the advanced solution of the Einstein equations the 
differential VLBI time delay in Eq. (29) is given by 

N 



Ait.M = 2^M.(l + K...(.^0) In ^^f , (34) 

^ V J r2a{s+2) + K-r2a{s+2) 

where s+i and s+2 are advanced times related to the times of arrival ti and t2 of radio signal 
to the first and second VLBI stations through the solutions of the advanced gravity cone 
Eqs. (25), (26) with e = +1. 

Taylor expansion of positions Xa{s+i) of the bodies around the coordinate times of arrival 
ti transforms Eq. (34) to 

= sf fl + ^) + ^ 2GM., (y.,xJC).(BxJC) 



Comparison of Eq. (35) with Eq. (28) elucidates that the signs of the second terms in their 
right sides are opposite. Hence, these terms originates due to the effect of propagation of 
gravity and can be used to measure its speed. This interpretation has been given by Kopeikin 
(2001) and is confirmed by calculations given in the present paper. Interpretations given by 
Asada (2002) and Will (2003) are erroneous and originates from misleading ideas. 
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7. Experimental Verification of the Retarded Nature of the Einstein Equations 

Let us recast Eqs. (28) and (35) to the following form (c.f. Eq. (12) from (Kopeikin 2001)) 



where the parameter S = — {e + 1) = c/ Cg — 1 with the speed of gravity Cg = c, if we operate 
with the retarded (e = —1) branch of the gravity null cone Eq. (1), and —Cg = c, if the 
advanced (e = +1) branch of the gravity null cone Eq. (1) is used. Hence, 5 = 0, if solutions 
of the Einstein gravitational field equations have a retarded nature, and 6 = —2, if the 
nature of the solutions of the Einstein equations is advanced. Furthermore, one could use 
a half-retarded plus half-advanced solutions of the Einstein equations in order to calculate 
the differential VLBI time delay. We have fulfilled such calculation and proved that it 
corresponds to the case of = oo in Eq. (36), that is 5 = — 1 (e = 0). Formally, it 
is equivalent to the case of the Newtonian theory in which the propagation of gravity is 
instantaneous with infinite speed (Kopeikin 2003). 

We have used the experimental data of the gravitational VLBI experiment on September 
8, 2002 (Fomalont & Kopeikin 2003) to distinguish between the three cases for the parameter 
6 = 0,-1,-2, and to measure the speed of gravity Cg. We found (Fomalont & Kopeikin 
2003) that 6 = -0.02 ± 0.19 and the speed of gravity Cg = (1.06 ± 0.21)c. Two cases 
6 = (red loop) and 6 = —1 (green loop) are shown in Fig 1. Our experiment convincingly 
demonstrates that only retarded solution 5 = is consistent with the observations proving 
that the speed of gravity Cg is equal to the speed of light c within the experimental errors. 




In 



Rla + KRia 
Rla + K-R2a 



(1+5) 



2GMj {vjxK)-{BxK) 
(Ru + KRu) 



(36) 
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Fig. 1. — While Jupiter moves in the sky from right to left the apparent position of the quasar 
J0842+1835 moves counterclockwise and completes a loop in the sky. If gravity propagates 
instantaneously, the apparent position of the quasar will move along the green circle 1-2-3-4- 
5 corresponding to the present positions 1, 2, 3, 4, 5 of Jupiter (white circles). The red curve 
l*-2*-3*-4*-5* is our theoretical prediction for the light deflection of the quasar based on 
general relativity and the retarded positions of Jupiter (brown circles). Measuring the size 
and shape of the observed loop (especially in the neighborhood of the point 3*) allows us 
to evaluate the difference between the speed of gravity and light and to answer the question 
whether the light is deflected by the present position of Jupiter (5 = —1) or by its retarded 
position (5 = 0) due to the flnite speed of gravity (Fomalont & Kopeikin 2003). 



